Manuscript prepared for Biogeosciences Discuss. 

with version 2.2 of the I5T|=X class copernicus_discussions.cls. 

Date: 8 November 201 1 


Filling-in of broad far-red solar lines by 
terrestrial fluorescence and atmospheric 
Raman scattering as detected by 
SCIAMACHY satellite measurements 

J. Joiner 1 , Y. Yoshida 2 , A. P. Vasilkov 2 , E. M. Middleton 1 , P. K. E. Campbell 3 , 

Y. Yoshida 4 , A. Kuze 5 , and L. A. Corp 6 

'NASA Goddard Space Flight Center, Greenbelt, MD, USA 

2 Science Systems and Applications, Inc., 10210 Greenbelt, Rd., Ste 400, Lanham, MD, USA 
3 University of Maryland, Baltimore County, Joint Center for Environmental Technology 
(UMBC-JCET), Baltimore, MD, USA 

4 National Institute for Environmental Studies (NIES), Tsukuba-City, Ibaraki, Japan 
5 Japan Aerospace Exploration Agency (JAXA), Tsukuba-City, Ibaraki, Japan 
6 Sigma Space Corp., Lanham, MD USA 

Correspondence to: J. Joiner (joanna.joiner@nasa.gov) 


1 



Abstract 


Global mapping of terrestrial vegetation fluorescence from space has recently been 
accomplished with high spectral resolution (v/Av > 35000 ) measurements from the 
Japanese Greenhouse gases Observing SATellite (GOSAT). These data are of interest 
because they can potentially provide global information on the functional status of veg- 
etation including light use efficiency and global primary productivity that can be used 
for global carbon cycle modeling. Quantifying the impact of fluorescence on the 0 2 -A 
band is important as this band is used for cloud- and aerosol-characterization for other 
trace-gas retrievals including C0 2 . Here, we demonstrate that fluorescence information 
can be derived from space using potentially lower-cost hyperspectral instrumentation, 
i.e., more than an order of magnitude less spectral resolution {u/Av ~ 1600 ) than 
GOSAT, with a relatively simple algorithm. As a demonstration, we use the filling-in of 
one of the few wide and deep solar Fraunhofer lines in the red and far-red chlorophyll- 
a fluorescence bands, the calcium II line near 866 nm, to retrieve fluorescence with 
the Scanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIA- 
MACHY) satellite instrument. Although the signal from vegetation fluorescence is ex- 
tremely weak at 866 nm, our results suggest that useful information may be obtained 
after adjustments are made to the observed spectra to correct for instrumental arti- 
facts. We compare fluorescence from SCIAMACHY with that retrieved at 758 and 
770 nm from similarly-corrected GOSAT data as well with the Enhanced Vegetation 
Index (EVI) from the MODerate-resolution Imaging Spectroradiometer (MODIS). We 
also show that filling-in occurs at 866 nm over barren areas, possibly originating from 
luminescent minerals in rock and soil. 


1 introduction 

Fluorescence from terrestrial vegetation has been measured in the laboratory and with 
ground- and aircraft-based instruments for several decades (see e.g., the review of 
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Meroni et al., 2009, and the many references therein). More recently, it has been 
demonstrated that fluorescence information can also be derived from satellite pas- 
sive sensors measuring backscattered sunlight at red and far-red wavelengths. These 
fluorescence measurements, particularly with global coverage from satellites, are of 
interest because the signal originates from the core complexes of the photosynthetic 
machinery where energy conversion of absorbed photosynthetically active radiation 
occurs. Chlorophyll fluorescence measurements therefore provide a means to esti- 
mate global instantaneous vegetation carbon-related processes. These processes are 
of major interest for precision farming, forest management, and assessment of the ter- 
restrial carbon budget including gross primary productivity (GPP) (e.g., Campbell et 
al., 2008; Damm et al., 2010; Frankenberg et al., 2011b). In addition, it is important 
to quantify the impact of fluorescence on atmospheric absorption bands, such as the 
O 2 -A band, that are used for cloud- and aerosol-characterization for other trace-gas 
retrievals including CO 2 (Frankenberg et al., 2011a). 

Fluorescence in vegetation occurs because a portion of the sunlight absorbed by 
chlorophyll cannot be used for carbon fixation. While most of this unused energy is 
given off as heat, a small fraction is re-emitted as fluorescence at longer wavelengths 
(lower energy) with respect to the excitation wavelength. The far-red and red fluores- 
cence generally occurs as a convolution of broad band emission with two peaks at 685 
and 740 nm, respectively, as shown in Fig. 1 (e.g., Meroni et al., 2009; Corp et al., 2003, 
2006). Most studies have shown that in high light conditions (i.e., in the afternoon when 
GOSAT measurements are made) and when plants are under stress, fluorescence is 
correlated with photosynthesis and light use efficiency (LUE) (e.g., Flexas et al., 2002; 
Louis et al., 2005; Meroni et al., 2008; van der Tol et al., 2009; Zarco-Tejada et al., 
2009; Amoros-Lopez et al., 2008). Variations in fluorescence and photosynthesis oc- 
cur when protective mechanisms intervene to prevent damage by harmful radicals that 
are formed under these conditions. 

Several studies have indicated that chlorophyll fluorescence supplies information 
content that is complementary to reflectance-based spectral vegetation indices includ- 
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ing the Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index 
(EVI), related to chlorophyll content, and the Photochemical Reflectance Index (PRI), 
related to changes in xanthophyll cycle pigments (Meroni and Colombo, 2006; Middle- 
ton et al., 2008; Rascher et al., 2009; Middleton et al., 2009). For example, Guanter 
et al. (2007) showed satellite observations with low variation of the NDVI within a field, 
while PRI and fluorescence ( F s ) showed significant variations. PRI and F s were not 
correlated, perhaps owing to the large dependence of PRI on the canopy structure. 
Zarco-Tejada et al. (2009) similarly showed no correlation between crown fluorescence 
and NDVI measured from airborne imagery. Meroni et al. (2008) showed that F s nor- 
malized by incoming photosynthetically active radiance (PAR) was correlated with LUE 
regardless of stress, while PRI showed a more complex relationship. Preliminary re- 
sults from Joiner et al. (2011a) showed that fluorescence derived from satellite mea- 
surements in some areas appears to lead EVI/NDVI by a month or so during green up 
and brown down periods. Both Joiner et al. (2011a) and Frankenberg et al. (2011b) 
showed areas and vegetation types for which fluorescence and vegetation indices dif- 
fered, particularly moderate to high latitude evergreen that maintains its pigment during 
winter inactivity. 

Space-based measurement of solar-induced chlorophyll fluorescence is challenging, 
because its signal (typically 1-5% red and far-red wavelengths) is small as compared 
with the much larger reflectance signal. Ground- and aircraft-based approaches have 
made use of the dark and spectrally-wide O 2 -A (—760 nm) and 0 2 -B (—690 nm) at- 
mospheric features to detect the weak fluorescence signal (see e.g., Meroni et al., 
2009). Deep solar Fraunhofer lines have also been used to detect fluorescence from 
vegetation using airborne and ground-based sensors. Initial studies (e.g., Plascyk and 
Gabriel, 1975; McFarlane et al., 1980; Watson, 1981; Sioris et al., 2003) focused on 
lines shortward of the red edge (— 700nm) where the terrestrial surface albedo is rela- 
tively low. More recently, Joiner et al. (201 1 a) and Frankenberg et al. (201 1 a,b) focused 
on longer-wavelength solar Fraunhofer lines that can be observed with space-based in- 
struments such as the currently operational GOSAT and the upcoming OCO-2. They 
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showed that fluorescence can be detected using Fraunhofer lines away from the chloro- 
phyll fluorescence peak even when the surface is relatively bright. 

Here, we build on that work making use of unique spectral measurements from 
the Scanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIA- 
MACHY) instrument aboard the European Space Agency (ESA) Environmental satel- 
lite (EnviSat). Spectral radiance measurements encompassing a wide solar Fraunhofer 
calcium (Ca) II line at 866 nm are used to retrieve an additive signal that fills in the line. 
These measurements are at the long wavelength tail of the chlorophyll-a fluorescence 
spectrum. 

Very few laboratory and ground-based measurements of vegetation fluorescence 
have been reported at wavelengths longer than 800 nm. Figure 1 displays measure- 
ments acquired using polychromatic excitation at wavelengths shorter than 665 nm 
with a so-called solar simulator after a five minute dark adaptation of the samples as 
described in Campbell et al. (2008). Chlorophyll fluorescence levels decrease signif- 
icantly from the dark-adapted maximum (F max ) to the steady state ( F s ) in 15-20 s, 
after which they remain relatively constant. Here, we show that F max and F s in the 
red and far-red peaks varies significantly, while similar and much lower values occur 
at the fluorescence spectrum tails (e.g., at 656 and 854-866 nm). Amoros-Lopez et 
al. (2008) similarly report spectral measurements from 650-850 nm showing a small 
amount of fluorescence (~0.1 mW/m 2 /nm/sr) at the longest wavelength when illumi- 
nated by sunlight through a filter at wavelengths shorter than about 600 nm (i.e., not 
full ambient intensity at leaf level). Other measurements (e.g., Chappelle and Williams, 
1987; Rosemaet al., 1991; Edner et al., 1994; Saito et al., 1998) suggest that the tail of 
the chlorophyll-a fluorescence bands and perhaps fluorescence from other compounds 
such as vitamin Ki (Chappelle et al., 1984) may produce a measurable signal at wave- 
lengths beyond 800 nm. Measurements have shown that the excitation efficiency of 
chlorophyll fluorescence is high in the blue and ultraviolet (UV) with a minimum around 
550 nm and another peak of high efficiency beyond about 650 nm (Chappelle and 
Williams, 1987; Rosema et al., 1991 ; Middleton et al., 2008). Therefore, solar-induced 
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fluorescence at long wavelengths may be higher than that measured by laser or in a 
partially-illuminated scenario. 

In addition to vegetation, minerals in rocks and soil may luminesce from blue to 
far-red wavelengths and contribute to the filling-in of Fraunhofer lines. Mineral detec- 
tion was a primary goal of early measurements using the Fraunhofer line discriminator 
method (e.g., Hemphill and Vickers, 1966; Watson, 1981; Stacy et at., 1984; Hemphill 
et al., 1988). For example, Watson (1981) found that luminescence at 489 nm from an 
airborne Fraunhofer line discriminator instrument had a good correlation with fluorite 
concentrations measured in soil samples collected along two traverses. 

A primary objective of this work is to explore the possibility of making fluorescence 
measurements from space with lower spectral resolution instrumentation than either 
the GOSAT interferometer or the OCO-2 spectrometer. If this can be demonstrated, 
then low cost instrumentation can potentially be used for such measurements with 
a relatively simple retrieval algorithm, i.e., one that does not require complex atmo- 
spheric correction. Of the wide and deep solar lines shown in Fig. 1, we focus on 
the 866 nm Ca II solar Fraunhofer line, because it is the only line with negligible con- 
tamination from atmospheric water vapor capable of detecting fluorescence at spectral 
resolutions of about 0.5 nm. However, this is clearly not an optimal wavelength region 
for measuring fluorescence owing to the weakness of the signal; instrumental artifacts 
can easily overwhelm the fluorescence signal. SCIAMACHY is the only satellite in- 
strument making hyperspectral measurements at 866 nm. Though not optimized for 
detecting fluorescence, it can be used for an exploratory study. 

The paper is organized as follows: Sect. 2 describes the satellite observations used 
here. Simulations of the Ca II line-filling at 866 nm from atmospheric rotational- and 
vibrational-Raman scattering as well as chlorophyll fluorescence are given in Sect. 3. 
The details of our retrieval approach, including methodology to correct for instrumental 
artifacts, are provided in Sect. 4. We show maps of monthly-mean fluorescence de- 
rived from 866 nm Ca II line-filling and a seasonal climatology derived from more than 
8 years of SCIAMACHY data in Sect. 5. In this section, we also compare and contrast 
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temporal and spatial variations of fluorescence (retrieved from both SCIAMACHY and 
GOSAT) and the MODIS enhanced vegetation index (EVI). Conclusions are given in 
Sect. 6. 

2 Satellite observations 
2.1 SCIAMACHY 

SCIAMACHY is a grating spectrometer that measures transmitted, reflected, and scat- 
tered sunlight in both limb- and nadir-viewing geometries from ultraviolet to near- 
infrared wavelengths (212-2386 nm) in eight separate channels. It was launched on 
the EnviSat platform in February 2002 into a sunsynchronous orbit with a descending 
node equator crossing time near 10:00. The instrument can make measurements in a 
number of different modes (e.g. Lichtenberg et al., 2006). In addition to Earth views in 
the limb and nadir, it can operate in a solar or lunar occultation mode. It also measures 
the solar irradiance. In this work, we primarily use channel 5 that covers wavelengths 
between 773 and 1063 nm at a spectral resolution of 0.54 nm. The signal-to-noise 
ratio for nadir measurements in this channel is between 1000 and 10000 (Gottwalk et 
al., 2006). We also examined spectra in channel 4 that covers 595-812 nm at 0.48 nm 
spectral resolution. 

We processed raw SCIAMACHY level IB data with the SciaLIc command-line tool 
software package (DLR, 2006). Using this package, we applied all available correc- 
tions and calibrations. Although the corrections do an excellent job of removing most 
artifacts that produce false filling-in of solar Fraunhofer lines, we found it necessary to 
further account for the remaining effects as discussed in Sect. 4. In App. A, we give a 
brief overview of the effects. More details can be found in Lichtenberg et al. (2006). 
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2.2 GOSAT 


GOSAT is a satellite mission designed to monitor the global distribution of the green- 
house gases C0 2 and CH 4 (Yokota et al., 2009). Jointly developed by the Japanese 
Ministry of the Environment (MOE), the National Institute for Environmental Stud- 
ies (NIES), and the Japanese Aerospace Exploration Agency (JAXA), GOSAT was 
launched on 23 January 2009 into a sun-synchronous orbit. It has a descending 
node equatorial crossing time near 13:00LT and a 3 day repeat cycle. GOSAT 
has two instrument packages: 1 ) The Thermal And Near-infrared Sensor for carbon 
Observation-Fourier Transform Spectrometer (TANSO-FTS) (Kuze et al., 2009) that 
measures backscattered solar radiation in three shortwave infrared (SWIR) regions, 
referred to as “bands”, centered at 0.76, 1 .6, and 2.0 pm with a nadir ground footprint 
of 10.5 km diameter and 2) The Cloud and Aerosol Imager (CAI) that contains four 
bands at 0.38, 0.67, 0.87, and 1.6 pm and footprints between 0.5 and 1.5km. 

Joiner et al. (2011a) and Frankenberg et al. (2011b) retrieved chlorophyll fluores- 
cence using TANSO-FTS band 1 . This band extends from approximately 758 to 775 nm 
with a resolving power (u/Av) of > 35000. The signal-to-noise ratio (SNR) for this band 
is > 300 for a typical scene radiance. 

There are also instrumental effects that can distort GOSAT spectra to produce a 
false filling-in of Fraunhofer lines. These include non-linearity from the analog-to-digital 
converter (ADC) (Kuze et al., 2011b) and the analog filter circuit in band 1 that can lead 
to effects such as a zero-level offset (Frankenberg et al., 201 1 b). In addition, calibration 
drift (Kuze et al., 2011a) may produce errors and drift in the absolute values of derived 
fluorescence. The latest available version (v130) of the level 1b data addresses some, 
but not all of these problems. Future versions will provide more corrections for these 
non-linearities. Only a limited amount of vl 30 data is currently available. Therefore, 
here we use a previous version of level 1 b data (v050) with an approach to account for 
instrumental artifacts that produce false filling in as described below. 
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2.3 MODIS Vegetation indices 


Below, we compare satellite-derived fluorescence with the Enhanced Vegetation In- 
dex (EVI) (Huete et at., 2002), a popular vegetation reflectance-based index that indi- 
cates relative greenness and is used to infer photosynthetic function. The EVI, derived 
from the MODIS instrument, can be found in the “Vegetation Indices 16-Day L3 Global 
1km” (MOD13A2) data set. We use data from the Aqua MODIS exclusively here. The 
Aqua satellite has an ascending node equator crossing near 13:30LT, similar to that 
of GOSAT. Although these data have been cloud-filtered and corrected for aerosol- 
corruption effects, we followed the methodology of Xu et al. (2011) to further remove 
residual atmospheric corruption in order to produce the best possible signals from veg- 
etation. 

3 Simulated filling-in at 866 nm 

Here, we simulate the effects of fluorescence as well as rotational- and vibration- 
Raman scattering on space-based observations of the Ca II line near 866 nm. We 
use the so-called KPNO2010 high spectral resolution solar irradiance reference spec- 
trum (Chance and Kurucz, 201 0) for these simulations and convolve the spectrum with 
triangular instrument line shape functions having different values of the full-width half- 
maximum (FWHM) including one similar to SCIAMACHY. 

3.1 Fluorescence 

Figure 2 shows simulated filling-in of the Ca II solar line at 866 nm for a spectrally con- 
stant fluorescence signal of 0.2 mW/m 2 /sr/nm for several different spectral resolutions 
including the limit at monochromatic resolution. Because Rayleigh optical thickness 
is low at these wavelengths, we did not include atmospheric scattering in these sim- 
ulations. The filling-in at SCIAMACHY resolution is small (<1%). However, modest 
increases in spectral resolution will increase the filling-in. 
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Figure 3 shows the surface-albedo dependence of the simulated filling-in of the Ca II 
line for the same amount of fluorescence (0.2 mW/m 2 /nm/sr) and spectral resolutions 
as in Fig. 2. The surface albedo dependence is significant at low surface albedos. 
However, the dependence flattens out at values typical of vegetated land (>-0.25). 

3.2 Rotational-Raman scattering 

We computed the filling-in of the Ca II line due to rotational-Raman scattering (RRS) 
using the scalar Linearized Discrete Ordinate Radiative Transfer code (LIDORT-RRS) 
of Spurr et al. (2008). Figure 4 shows the filling-in of the Ca II line as in Fig. 2 but 
from rotational-Raman scattering (RRS) only. Although the amount of atmospheric 
scattering is relatively small at these wavelengths, the RRS filling-in is apparent owing 
to the large depth of this solar line as similarly shown by Joiner et al. (2011a) for the 
K I line at 770 nm. The spectral signature is very similar to that of line filling due to 
fluorescence. The filling-in due to RRS is about a factor of 6 less than that due to 
a fluorescence of 0.2 mW/m 2 /nm/sr at SCIAMACHY resolution. This factor increases 
somewhat at higher spectral resolutions. 

Figure 5 shows that the RRS filling-in dependence on surface albedo is relatively flat 
for albedos >-0.3. There is a small but significant dependence on both surface pres- 
sure and solar zenith angle (SZA). Figure 6 shows the SZA dependence in more detail 
as well as the dependence on satellite view zenith angle (VZA) and surface pressure. 
The SZA dependence is relatively small for SZA < 60°, but increases rapidly for greater 
angles. The VZA dependence is small for typical angles viewed by SCIAMACHY. 

3.3 Vibrational-Raman scattering 

Vibrational Raman scattering (VRS) is commonly neglected in UV and visible retrieval 
algorithms as it is small compared with rotational-Raman scattering (RRS) (Burrows et 
al., 1996). However, the fractional contribution of VRS is larger in the near-InfraRed 
(NIR) owing to higher values of solar radiation at excitation wavelengths as compared 
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with emission wavelengths. The VRS frequency shifts between excitation and emis- 
sion are 2331 and 1555 cm -1 for N 2 and O 2 , respectively, assuming pure vibra- 
tional transitions of the molecules (Heise and Schrotter, 1995). We assessed the VRS 
contribution to the filling-in of the Ca II Fraunhofer line using the single scattering ap- 
proximation of the radiative transfer equation. Along with photons only scattered in the 
atmosphere, we also accounted for photons first reflected from the surface and then 
scattered in the atmosphere. 

Absolute differential Raman cross sections of the vibrational line at 488 nm at 90° 
measured by different authors for N 2 vary from 3.3 x icr 31 to 5.6 x 10~ 31 cm 2 /str-mol 
(Fenner et al., 1973). Here, we use a value of 5.5 x to -31 cm 2 /str-mol that is consistent 
with data from Heise and Schrotter (1995). The ratio of O 2 to N 2 cross sections at 488 
nm is assumed to be 1.3 (Fenner et al., 1973). We adopt a general form of phase 
functions for both the unshifted Cabannes line and wavelength shifted VRS lines as 
suggested by Danichkin et al. (1982). For the central Cabannes line in air, the depo- 
larization ratio for unpolarized light is assumed equal to 0.00734 (Chance and Spurr, 
1997). For VRS lines, we use values recommended by Danichkin et al. (1982), i.e., 
0.043 and 0.099 for N 2 and O 2 , respectively. 

Figure 7 shows the VRS filling-in of the 866nm Ca II Fraunhofer line as a function of 
surface albedo. The maximum value of about 0.45% occurs for a non-reflecting surface. 
VRS filling-in rapidly decreases with increasing surface albedo. For typical values of 
surface albedo over land (>0.2), the filling-in only weakly depends on surface albedo 
and does not exceed 0.02% for SZA< 60°. Therefore, we neglect VRS in the remainder 
of this study. 
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4 Retrieval methodology 

4.1 Fitting windows 

The use of solar Fraunhofer lines substantially simplifies a space-based retrieval of F s 
as compared with the 0 2 -A band as described by Joiner et al. (2011a) and Franken- 
berg et al. (2011a); atmospheric correction is unnecessary as long as observations 
are made at wavelengths not affected by atmospheric absorption. Joiner et al. (2011a) 
used the potassium (K) I line near 770 nm to derive fluorescence using GOSAT’s 
TANSO-FTS. This line falls between several weak 0 2 lines and can be observed in rela- 
tive isolation with GOSAT’s high spectral resolution. Here, we derive fluorescence from 
GOSAT with the same fitting window used by Joiner et al. (201 1 a) (769.90-770.25 nm). 
We use a second fitting window between 758.45 and 758.85 nm and similar to Franken- 
berg et al. (2011b) scale the results from that window and add them to the other to 
increase the signal-to-noise ratio. 

Of the solar lines shown in Fig. 1 (656, 850, 854, 866 nm), only the 866 nm Ca 
II line is relatively free of contamination by atmospheric water vapor absorption. We 
attempted fluorescence retrievals using all of these lines including the Ha line near 656 
nm in SCIAMACHY channel 4. We found that it was extremely difficult to disentangle 
the fluorescence filling-in from water vapor absorption. Therefore, in the remainder of 
this work, we focus solely on deriving fluorescence using the 866 nm Ca II solar line 
from SCIAMACHY with the spectral band 863.5-868.5 nm. 

4.2 General approach 

We use a variation on the approach implemented by Joiner et al. (2011a) to jointly re- 
trieve a reflectivity-related parameter K and an additive signal F s presumably produced 
by some type of fluorescence. We use the above-mentioned spectra fitting windows 
with the following simplified model for the observed Earth spectral radiance 1(A) that 
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assumes negligible atmospheric absorption and scattering: 

I(\) = KE(\)*+F S , (1) 

(Joiner et at., 2011a) where 23(A)* is a reference spectrum containing no fluorescence. 

As in Joiner et al. (2011a), we assume a constant wavelength dependence of F s 
over our relatively small fitting windows. We fit both constant and linear in wavelength 
coefficients for K (i.e., K = K 0 + K\ [A - A 0 ] ) for the larger SCIAMACHY fitting win- 
dow, but only a constant for the smaller GOSAT window. The spectral dependence of 
atmospheric, cloud, and aerosol scattering is essentially constant over the relatively 
small narrow fitting windows used, so the neglect of it will have an insignificant effect 
on the retrieval of fluorescence. Similarly, calibration of the reference spectrum 23(A)* 
does not affect the derived fluorescence, since any offset between the reference and 
observed spectra is absorbed in the parameter K. An error in absolute calibration of 
1(A), however, will produce an error in the retrieved F s . 

Similar to the approach of Joiner et al. (2011a), we use a standard weighted least 
squares fitting procedure. All wavelengths are weighted equally except for a few wave- 
lengths within the GOSAT fitting window that are not used because they are affected 
by O 2 absorption. 

4.3 Use of reference spectra 

The main difference between our approach and that of Joiner et al. (2011a) is that 
we use spectral radiance measurements made over the cloudy ocean as a reference 
rather than a measured or computed solar irradiance spectrum. No further corrections 
to observed spectra are made. Our retrieval approach is insensitive to errors in the 
instrument line shape function or reference solar spectra as these are not used. 

The main criterion for reference spectra used to derive terrestrial fluorescence is that 
they should represent the spectral shape of a solar Fraunhofer line without terrestrial 
fluorescence. There are several advantages of using cloudy Earth radiance spectra as 
a reference to derive the terrestrial fluorescence signal as opposed to solar irradiance 
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spectra (measured or computed). Firstly, the use of Earth spectra minimizes the so- 
called undersampling problem and eliminates errors due to inaccuracy of the assumed 
instrument line shape function. The undersampling problem and a proposed correction 
method have been discussed at length by Chance et al. (2005). However, we found 
that the proposed correction method did not achieve the desired accuracy owing to 
uncertainties in the instrument line shape function. Because the wavelength shifts are 
very small between ocean-based reference spectra and the land-spectra of interest, 
problems due to undersampling are essentially eliminated. Corrections used by Joiner 
et al. (201 1 a) to account for undersampling when referencing to observed solar spectra 
are not necessary. 

Secondly, referencing to ocean spectra with similar continuum radiance values al- 
lows us to significantly reduce the errors from straylight, dark current, zero-level off- 
set, and other instrumental effects described above that can produce a false filling-in. 
These effects should be present in the ocean reference spectra at approximately the 
same level. Therefore, the effects should be accounted for when deriving fluorescence 
using Eq. 1. 

Finally, using ocean spectra at similar radiance levels and SZAs as a reference will 
also implicitly account for some of the filling-in due to RRS. RRS depends upon sur- 
face albedo as well as the mean pressure of the scattering surface. Here, we are using 
cloudy data at similar SZAs and continuum radiance values. Studies of the so-called 
optical centroid pressure (approximately, the reflectance-weighted mean pressure) of 
clouds over ocean at moderate reflectivities typical of vegetated land show that the 
clouds occur mainly between about 775 and 875 hPa (Joiner et al., 2011b). Deep 
convective and frontal clouds that have lower optical centroid pressures typically oc- 
cur only at higher Lambertian-equivalent reflectivities (>~ 0.6). Therefore, the main 
difference in RRS filling-in between the cloudy ocean and land spectra is due to the 
different effective surface pressures. As shown in Fig. 5, this difference is small, but 
not negligible for the expected pressure differences. We must bear this difference in 
mind when examining our results. For regions with surface pressures higher than about 
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900 hPa, there may be unaccounted for filling-in due to RRS. However, for high-altitude 
land regions (< 775hPa) the reference ocean spectra may overestimate RRS filling-in. 

Frankenberg et al. (2011b) used a different approach to account for instrumental arti- 
facts; they utilized measurements over Antarctica to develop a parameterization for the 
false fluorescence produced by the zero-level offset in GOSAT data. They showed gen- 
erally increasing amounts of false fluorescence with radiance levels with some wiggles 
and a downturn at the highest radiances in the S polarization. We found similar features 
using cloudy ocean data. With our approach, we are able to cover the full range of ob- 
served radiance values and see a downturn in the false fluorescence or filling-in at the 
highest radiance levels also in the P polarization. Our utilization of cloudy ocean data 
does not necessitate the use of high SZA data to cover the full range of observed con- 
tinuum radiance values; filling-in from RRS increases rapidly at high SZAs as shown 
above, although the imprint of RRS is not obvious in the results of Frankenberg et al. 
(2011b). 

Frankenberg et al. (201 1 b) discarded GOSAT P polarization data because the errors 
in that polarization were time-dependent. We likewise find such time-dependent errors 
in both GOSAT and SCIAMACHY data. Because cloudy ocean data are available year 
round, we are able to account for these errors and use all data including GOSAT P 
polarization, thereby significantly reducing the effects of instrumental noise. 

In addition, cloudy ocean data span all latitudes. Therefore, we are able to account 
for instrumental effects that vary systematically with orbital phase; our approach is 
able to detect and correct for subtle but significant latitudinally-dependent artifacts in 
both GOSAT and SCIAMACHY data. We find that the latitudinal-dependence of these 
artifacts also varies with season. 

Here, we binned potential reference spectra according to several criteria discussed 
below. For each bin, we averaged multiple spectra to produce a composite reference 
spectrum. To reduce any filling-in signal from the ocean, e.g., from ocean Raman 
scattering (Vasilkov et al., 2002), RRS, VRS, or fluorescence, we only use pixels with 
reflectances > 0.2 for SCIAMACHY. For GOSAT we used only pixels with a cloud frac- 
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tion (derived from the CAI) of >20%. 

4.3.1 SCIAMACHY binning 

For SCIAMACHY, we binned spectra separately for each of the two primary pixel expo- 
sure times and for each scan position within a nadir block according to the continuum 
radiance value and the wavelength shift (with respect to a fixed solar spectrum). We 
found that the observed filling-in over ocean sometimes changes abruptly and that it 
also drifts slowly over time. We therefore compute a set of reference spectra for each 
day of observations. If the standard deviation of the filling-in over ocean for a given day 
exceeds a threshold value, we eliminate that day from our sample. 

We use a bin size of 0.1 for (normalized) radiance, represented by .Rcos(SZA) where 
R is the reflectance. We found that the wavelength shifts varied systematically with 
latitude, so binning by both wavelength shift and reflectivity implicitly accounts for filling- 
in features such as RRS that are function of solar zenith angle as well as instrumental 
artifacts that vary systematically with orbital phase. Because the wavelength shifts vary 
with time, we compute the daily mean and standard deviation of wavelength shift and 
use 8 bins with intervals of half the standard deviation of the wavelength shift (typically 
-0.007 nm). 

4.3.2 GOSAT binning and additional corrections 

For GOSAT, we used bin sizes of 0.05 for normalized radiance, 20° for latitude, and 
0.001 nm for wavelength shift. We use only data processed with the H gain range as 
there are not enough data from the M range to perform a reliable correction. Therefore, 
data over the Sahara and in parts of central Australia have not been processed. For 
GOSAT, we retrieve F a separately for the S and P polarizations (separate corrections 
applied) and combine them. We also apply a correction factor to the radiances to 
account for instrument degradation as suggested by Kuze et al. (2011a). 
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4.4 Cloud filtering 


Because clouds do not alter the spectral structure used to retrieve F s , we can tolerate 
small amounts of cloud within our satellite footprints; the main effect of clouds will be 
to shield a portion of the F s signal from satellite view. We remove significantly cloud- 
contaminated data from our samples by two different methods. For SCIAMACHY, we 
use the concept of the cloud radiance fraction (CRF). CRF is defined as the fraction of 
pixel radiance scattered by cloud particles. The CRF is a convenient quantity for cloudy 
pixel screening, because it can be computed directly for large pixels and collocation 
with a subpixel imager is not necessary. The CRF is given by 

CRF = /^, (2) 

^obs 

where / is the geometric cloud fraction, and I obs and Im are the observed and bare- 
cloud (no surface reflection) radiances, respectively. Here we use the so-called mixed 
Lambertian-Equivalent Reflectivity model to estimate the CRF. In this model, I M is 
specified as a Lambertian surface with an equivalent reflectivity of 0.8, and / is set 
equal to an effective cloud fraction / eff that is derived using the relationship 

^obs = /eff-^cld "F (1 /eff)Tclr- (3) 

(e.g., see Stammes et al., 2008, for an overview of the approach that has been used to 
retrieve the cloud optical centroid pressure at NIR through UV wavelengths). In short, 
the selection of 0.8 as the cloud reflectivity produces an effective cloud fraction that 
is generally less than the geometric cloud fraction. Absorption and/or scattering from 
beneath a thin cloud is accounted for by the second term in Eq. 3. Here, we neglect 
atmospheric scattering as it is very small at 866 nm. To compute I clr , we use the black- 
sky 16-day gridded filled-land surface albedo product from Aqua MODIS (MOD43B3) 
at 865 nm (Lucht et al., 2000). We include only pixels with CRF <10% in our sample. 
High surface albedos over vegetated land and uncertainty in the SCIAMACHY absolute 
calibration relative to that of MODIS cause this check to allow a large fraction of the 
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data into our sample, removing only highly clouded pixels. We did not find a significant 
dependence of F s on the CRF and more stringent limits on cloud contamination did not 
significantly alter the derived spatial patterns of F s . 

Owing to uncertainty in the absolute calibration of GOSAT, here we use the cloud 
fraction derived from the GOSAT CAI. As in Joiner et al. (2011a), we removed pixels 
with estimated cloud fractions >10%. Joiner et al. (2011a) found that this provided 
similar results to those obtained when the CRF from the Ozone Monitoring Instrument 
(OMI) onboard the National Aeronautics and Space Administration (NASA) Aura satel- 
lite was used for GOSAT cloud detection. Results did not depend substantially on the 
exact values of the thresholds used here for both GOSAT and SCIAMACY. 

4.5 Quality control 

To create gridded data sets, we tried different filtering schemes for quality control. 
For SCIAMACHY, we included all data passing gross quality assurance checks on the 
filling-in, radiance residuals, and wavelength shifts. We tested whether eliminating the 
initial scan positions within a nadir-viewing block of observations affected the monthly 
means. This did not significantly affect the monthly means so all swath positions are 
included here. 

For GOSAT, we tested various filtering methods such as utilitizing consistency be- 
tween S and P polarizations and radiance residuals. In general, these checks ap- 
peared to remove noisy, but good data that was generally helpful in reducing random 
errors. Therefore, we did not apply any filtering to GOSAT data. 

4.6 Scaled fluorescence 

When comparing fluorescence with vegetation indices, we use a quantity called 
“scaled-F” (Joiner et al., 2011a), defined as the retrieved fluorescence F s divided by 
cos(SZA). This scaling roughly accounts for variations in F s due to the incoming (clear- 
sky) PAR. 
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4.7 Uncertainties 


Uncertainties for GOSAT retrievals were calculated by Joiner et al. (2011a) assuming 
only random errors (i.e., instrumental noise). However, errors from unaccounted-for 
systematic instrumental errors were quite significant. The absolute GOSAT F s values 
in this work are smaller (by approximately a factor of 2) than those derived by Joiner et 
al. (2011a), owing to the new corrections applied. 

Because unaccounted-for systematic errors in SCIAMACHY data could be signifi- 
cant and their magnitudes are unknown, it is difficult to quantify the uncertainties in the 
derived scaled-F. Many SCIAMACHY observations are available for averaging as com- 
pared with GOSAT which helps reduce random errors. However, systematic errors are 
expected to remain. For example, in addition to residual errors left over from imperfect 
corrections, the South Atlantic anomaly (SAA) has been known to affect measure- 
ments over South America. Here, we provide standard errors and standard deviations 
for gridbox averages to indicate subgridbox variability. We note that there may still be 
significant errors in the derived absolute values of F s and scaled-F. 

5 Results and Discussion 

5.1 Spatial variations in fluorescence 

Figure 8 shows retrieved gridded monthly mean scaled-F for July and December 2009 
derived from GOSAT at 770 nm and 758 nm (0.696*758 nm + 770 nm) and from SCIA- 
MACHY at 866 nm. The values from SCIAMACHY at 866 nm are significantly smaller 
than those at 770 nm, as expected. The expected seasonal variation in scaled-F is 
seen by both sensors, namely, higher Northern Hemisphere terrestrial activity in July 
versus higher activity in the Southern Hemisphere in December. Like the EVI, the re- 
trieved scaled-F is sensitive to the amount of green biomass contained within the sen- 
sor field-of-view or fractional amount of intercepted PAR. Similar patterns are shown in 
scaled-F by both SCIAMACHY and GOSAT, such as low values over central Brazil in 
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July and over the northeast of Brazil in December. 

Subtle differences in spatial variations are shown between GOSAT- and 
SCIAMACHY-derived scaled-F. SCIAMACHY shows somewhat less fluorescence at 
high northern latitudes in the boreal summer as compared with tropical regions. There 
could be several reasons for this. Firstly, there may be remaining errors due to instru- 
mental effects. Secondly, SCIAMACHY and GOSAT observations are made at different 
local times; SCIAMACHY observations are made near 10:00 while those of GOSAT are 
made about 2.5 hours later, just after noon. The shapes of the fluorescence spectra 
may differ for different types of vegetation as well as different magnitudes and spectral 
shapes of the incoming solar light. For example, more UV light is absorbed in rela- 
tively clear skies at the surface in the tropics where total column ozone and SZAs are 
low. Plants at these latitudes have more photo-protective mechanisms in place that 
may produce different fluorescent spectral responses. Compounds within the vegeta- 
tion besides chlorophyll-a may also contribute to the observed fluorescence. Finally, 
canopy structure may influence the fluorescence spectra. Multiple scattering within a 
deep canopy may enhance fluorescence, producing a different spectral response than 
a sparse canopy (Rosema et al., 1991). 

Figures 9-10 show a monthly climatology of fluorescence derived from over 8 years 
of SCIAMACHY data (January 2003-July 201 1 ) gridded at 0.5° x 0.5° resolution and a 
similarly gridded climatology of EVI from Aqua MODIS for comparison. An animated 
version of the climatology is available in the supplemental material. The general pat- 
terns from SCIAMACHY are consistent with fluorescence related to vegetation. Inter- 
esting differences between scaled-F and EVI are shown clearly here as well as in Fig. 8 
including moderate values of EVI over parts of the northeast and southeast US, Eu- 
rope, and southeast Asia year in the late autumn through early spring when scaled-F 
is close to zero. The values of scaled-F are more distinctly lower over parts of east and 
northeast Brazil from May through November as compared with the EVI. Scaled-F also 
falls within the moderate range for larger areas of Australia than EVI during its sum- 
mer months (November-February). These differences are further explored in the next 
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subsection. Note that we show scaled-F here for solar zenith angles up to 80°. In the 
winter hemisphere at the highest latitudes shown, scaled-F is likely to be overestimated 
owing to RRS that is not fully accounted for. 

5.2 Seasonal variations in fluorescence and EVI 

In Fig. 12, we examine the climatological seasonal cycle of fluorescence over several 
geographic regions shown in Fig. 1 1 . Within each region of interest, we further strat- 
ify the data based on the dominant International Geosphere Biosphere Programme 
(IGBP) land cover class. We show examples for several different biome types, using 
the same scales in all panels to highlight differences between scaled-F and EVI for the 
different regions and vegetation types. 

Over the pine forests in the southeast US (Fig. 12, panel 1), there is little lag between 
fluorescence and EVI. Deciduous and mixed forests over the southeast US showed pat- 
terns similar to those of evergreen needleleaf with little lag of the EVI (not shown). For 
this and for many other moderate to high latitude regions, values of the EVI reach min- 
ima between 0.2 and 0.3 during winter months, while scale-F is essentially zero. The 
one exception is southwestern Australia croplands (panel 8), where the EVI reaches a 
minimum value near 0.1 . 

Closed shrubland over northeast Brazil (not shown) exhibited a seasonal cycle very 
similar to open shrubland (panel 3), though with higher values of both scaled-F and 
EVI. In the open shrubland, both scaled-F and EVI are low in August through November 
following the wet season green-up that shows peaks in fluorescence and EVI in April. 
This contrasts sharply with the more heavily forested areas of the central Amazon 
(panel 2) at similar latitudes where both the EVI and scaled-F remain relatively high 
year round. 

Scaled-F appears to lead the EVI during the green up periods over croplands in the 
boxed areas of India and Thailand similar to the initial GOSAT results of Joiner et al. 
(2011a). There is a slightly earlier peak in scaled-F for near-equatorial Thailand (June) 
as compared with a later peak in India (July) presumably owing to latitude differences. 
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The lag in the EVI is also much greater during greenup over India (panel 4) than over 
Thailand (panel 5). In contrast, there is little seasonal variability in the EVI over forested 
regions of equatorial Indonesia (panel 6), while the fluorescence shows more variability 
with a minimum around January. 

In the forested parts of eastern Australia (panel 7), we see a clear seasonal cycle 
in scaled-F with a peak in the months December through February. The EVI shows 
a weaker seasonal variation. In contrast, the EVI shows a sharp peak in August and 
September over the croplands of southwest Australia (panel 8) while scaled-F shows 
much less variation. 

Higher northern latitudes typically show more seasonal variation and relatively high 
magnitudes of EVI in summer as compared with the tropics, while the SCIAMACHY 
scaled-F results suggest show lower peak magnitudes at these latitudes as compared 
with the tropics. An example is shown over areas of mixed forests and grass in France 
(panel 9). It is not known whether the relatively low scaled-F values obtained in the 
norther latitude summer are the result of remaining instrumental bias or represent a 
biophysical feature. 

The savannas of Africa (panels 1 0,1 1 ) show distinct seasonal cycles in both scaled-F 
and EVI. Again, scaled-F appears to lead the EVI, during both green-up and brown- 
down periods. The amount of lag varies with region. 

Finally, we show an example of deciduous broadleaf in the southern hemisphere 
(Bolivia, panel 12). Here again we see strong seasonal cycles in both scaled-F and 
EVI with a slight lag in the EVI during green-up. This region contains a wide variety 
of biome types. In this region, we also examined areas primarily consisting of open 
shrubland and grasslands (not shown). In these areas, we found very low (near zero) 
values of scaled-F as well as low values of EVI (~0. 1-0.2), sharply contrasting with 
the areas primarily covered by deciduous and evergreen (also not shown, but similar 
to deciduous) broadleaf vegetation. 

The apparent lags of EVI behind the fluorescence signal in several of our examples 
may have an alternative interpretation related to the fluorescence quantum yield that 
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is affected by environmentally-induced stress and other factors. The combination of 
F s and EVI (or the related parameter fraction absorbed PAR or FPAR) may enable 
partitioning of ecosystems for the purpose of improving monthly global GPP estimates. 
An additional measurement of F s in the 680 nm fluorescence feature may further aid 
in this partitioning. 

5.3 Unexplained filling-in (not due to vegetation fluorescence) 

Figure 13 shows a similarly derived monthly climatology of reflectivity at 866 nm. Taken 
in conjunction with Figs. 9 and 10, this indicates that fluorescence is loosely correlated 
with reflectivity in vegetated areas as would be expected if the filling-in of the Ca II 
line has a biogenic source. However, there are also areas of high reflectance with low 
fluorescence such as over the Sahara and parts of the Saudi Arabian peninsula. This 
demonstrates that the filling-in is not entirely related to reflectivity as might be the case 
if the signal was entirely due to an instrumental artifact. 

We note a significant filling-in (retrieved as fluorescence) over parts of the Sahara 
desert (most apparent in May-September) and the Saudi Arabian peninsula where 
vegetation is sparse and a significant signal from chlorophyll is not expected. These 
signals sometimes, but not always, occur over very high albedo surfaces. Examination 
of residual filling-in over Greenland shows residual filling-in over only the brightest area 
that may not have been well characterized by the cloudy spectra. Cloudy spectra 
with high radiance values occur primarily in deep and frontal convection where cloud 
pressures are low. This may lead to an underestimation of filling-in from RRS when 
these spectra are used as a reference. 

Filling-in over barren regions may also be produced by luminescent minerals in soil 
and/or rock. There are many types of minerals with impurities that luminesce when 
excited by UV and visible light (Gaft et al., 2005). For example, samples of fluorite 
(CaF 2 ) containing Nd 3+ (Neodymium) centers show laser-induced luminescence with 
peaks at 866 and 795 nm (Gaft et al., 2005). Similarly, Dy 3+ (Dysprosium) centers give 
narrow peaks at 765 and 778 nm. 
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High filling-in/fluorescence values are also frequently seen around continental and 
island coastlines, particularly on the Eastern edges. No such features are present 
in the reflectance climatology shown in Fig. 13. These high values may be due to 
a spatial aliasing effect or some other instrumental artifact such as distortion of the 
spectral response function in the presence of large radiance gradients within a satellite 
footprint. However, we observed neither high radiance residuals nor large wavelength 
shifts that may accompany such an occurrence. 

6 Conclusions and ongoing work 

Our results indicate that terrestrial fluorescence filling-in of the 866 nm Ca II line can be 
detected using hyperspectral instruments (spectral resolutions of the order of tenths of 
a nm) such as SCIAMACHY. Distinct and similar seasonal cycles of fluorescence are 
retrieved with SCIAMACHY at 866 nm and GOSAT at 770 nm. This finding is important 
because it suggests that fluorescence information may be retrieved using potentially 
low-cost air- and space-borne instrumentation with a relatively simple algorithm that 
does not rely upon complex atmospheric correction. The detected filling-in signals pre- 
sumably originate primarily from chlorophyll-a in vegetation though other compounds 
in vegetation may contribute. In addition to rotational-Raman scattering, luminescent 
minerals in soil and/or rock may also contribute to red and far-red solar Fraunhofer line 
filling. 

Owing to the large number of observations collected over more than 8 years of op- 
eration, SCIAMACHY data can be averaged to provide climatological spatial patterns 
of fluorescence at a higher resolution than is possible with the current set of GOSAT 
observations, even though the intrinsic pixel spatial resolution of SCIAMACHY is lower 
than that of GOSAT. The derived spatial patterns and seasonal cycles in fluorescence 
are similar to those of the MODIS EVI. However, important differences between scaled- 
F and EVI are noted. For example, fluorescence appears to lead EVI during greenup 
and brown down periods in some regions and for some vegetation types similar to re- 
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suits obtained previously by Joiner et al. (2011a) from GOSAT. A more distinct shut-off 
of activity in winter and dry months is shown in fluorescence in some areas as com- 
pared with EVI. 

More work needs to be done to further refine fluorescence retrievals from satellite 
instruments that were not designed to measure small Fraunhofer line filling-in signals. 
Owing to the extremely small fluorescence signal measured with instruments that have 
relatively large systematic errors, our reported absolute values and temporal/spatial 
variations may still contain residual errors. However, the approach we developed ap- 
pears to be capable of detecting and correcting for complex, time-varying instrumental 
artifacts. 

The influences of canopy structure and other large-scale parameters, such as chloro- 
phyll content, need to be further investigated in order to understand their effect on 
large-pixel satellite-retrieved fluorescence (e.g. Rosema et al., 1991 ; Louis et al., 2005; 
Meroni and Colombo, 2006; Middleton et al., 2008; van der Tol et al., 2009; Zarco- 
Tejada et al., 2009; Rascher et al., 2009; Damm et al., 2010). Indeed, it has been 
suggested that when vegetation is not under stress and the variations in fluorescence 
yield are small, fluorescence flux from the 740nm emission feature, that is not affected 
by reabsorption, contains direct information about absorbed PAR (Moya et al., 2004). 
In addition, the ratio of fluorescence from the 740 and 680 nm features expresses 
information about chlorophyll concentration (e.g., Gitelson et al., 1999; Saito et al., 
1998; Louis et al., 2005). Therefore, additional measurements of fluorescence at other 
wavelengths will be useful for disentangling the yield information from the chlorophyll 
content. In order to more fully exploit satellite-derived fluorescence measurements 
for deriving relationships with key carbon-related parameters such as gross primary 
production (GPP) or light-use efficiency (LUE), they should be combined with other ob- 
servable quantities such as the reflectance-based vegetation indices including the PRI. 
Additional laboratory and field measurements of fluorescence and reflectance spectra 
for a variety of vegetation types, canopy structures, illumination, and environmental 
conditions will help to improve the interpretation and use of satellite-derived fluores- 
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cence information. 

The upcoming launch of the Orbiting Carbon Observatory-2 (OCO-2) (Crisp et al., 
2004) will facilitate the use of solar Fraunhofer lines for fluorescence retrievals at sim- 
ilar wavelengths near the O 2 -A band (Frankenberg et al., 2011a) and provide more 
opportunities for fluorescence measurements. The FLuorescence Explorer (FLEX) 
(Rascher, 2007; European Space Agency, 2008), an ESA Explorer Mission selected 
for Phase A/BI in early 2011, plans to utilize the O 2 -A and B bands for chlorophyll 
fluorescence retrievals (Guanter et al., 2010). FLEX would provide fluorescence in- 
formation at a higher spatial resolution than current satellite sensors that are able to 
measure fluorescence but were not designed and optimized for that purpose. 
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Appendix A Instrumental effects in SCIAMACHY data that may produce false 
filling-in 

Al Memory effect 

In channels 1-5, it was found that the signal deviates from a linear response. The 
deviation depends upon the signal level of the previous readout. It also depends on the 
signal level including analog offset and dark current. Since there is a dependence on 
the previous signal for a large range of detector fillings, artificial spectral features can 
be introduced into the measurements. 
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We found systematic biases in Fraunhofer line filling as a function of position within 
a block of nadir measurements. The first several positions (following the limb-viewing 
mode) have the largest biases. Because of repeating orbital patterns, if this effect is 
not accounted for, a cross-track striping will be present in the observed filling-in, even 
after averaging a month of data. 

A2 Dark current 

There are two components that contribute to the dark signal in channels 1 -5: the analog 
offset and the leakage current. The former does not depend on integration time while 
the latter does. We found different biases in the filling-in for the two primary pixel 
exposure times (PETs) used for nadir measurements in channels 4 and 5 that may be 
a result of unaccounted for dark signals. 

A3 Wavelength calibration and undersampling 

Due to the Doppler and other effects, there are small wavelength shifts in Earth spectra 
with respect to the measured solar spectra. We found that these wavelength shifts vary 
systematically with latitude. Wavelength shifts between observed and reference spec- 
tra must be accurately accounted for when deriving a relatively small filling-in signal. 

A4 Straylight 

Stray light may be classified as two types: spectral and spatial. Both of these can 
distort observed spectra. Spatial stray light mainly affects the limb measurements of 
SCIAMACHY, but may also produce errors in a fluorescence retrieval in the nadir mode. 
Spectral stray light may be caused by diffuse reflection that adds signal to all detector 
pixels in a given channel, thus mimicking fluorescence. The largest expected uniform 
stray light fraction in SCIAMACHY is 0.14% (Lichtenberg et at., 2006). 

So-called “ghost” stray light is caused by reflection from one wavelength to another. 
It does not add signal to all pixels equally, but can distort the spectrum. The correction 
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method is accurate to around 25%, leaving at most 1% stray light in the spectrum. 

A5 Calibration 


Various issues with SCIAMACHY’s radiometric calibration are discussed by Lichten- 
berg et al. (2006). These may affect the absolute values of derived fluorescence. In 
addition to radiometry, the applied calibration scheme may introduce spectral features 
and thus produce false filling-in. 
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Fig. 1. Measured fluorescence spectra of excised corn leaves (Zea Mays L.) in the laboratory, 
with simulated solar illumination at full mid-day sun intensity as measured below a long-pass 
filter blocking wavelengths longer than 665 nm, showing maximum fluorescence (F max ) after 5 
min. of dark adaptation, steady state fluorescence (F s , achieved in less than 0.5 minute after 
illumination), and a sample solar spectrum from Thuillier et al. (2003, 2004) (1 nm spectral 
resolution). 
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Fig. 2. Simulated spectral filling-in around the Ca II solar line from fluores- 
cence^. 2 mW/m 2 /sr/nm, SZA=40°, nadir, and surface albedo=0.3 (RRS not included) for dif- 
ferent instrument resolutions (triangular slit function) with full-width half-maximum (FWHM) as 
indicated including that of SCIAMACHY (SC1A). 
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Fig. 3. Filling-in due to fluorescence=0.2 mW/m 2 /nm/sr, SZA=20°, and nadir as a function of 
surface albedo for different instrument spectral resolutions. 
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Fig. 4. Simulated filling-in spectra due to rotational-Raman scattering (RRS) only at SZA=40 
nadir, and surface albedo=0.3 for different instrument spectra resolutions. 





Surface albedo 


Fig. 5. Filling-in of Ca II line at 866 nm due to rotational-Raman scattering as a function 
of surface albedo for two solar zenith angles (SZA) and two surface pressures (P s ) showing 
simulated difference between low clouds over ocean (800 hPa) and a low lying land area (1000 
hPa). 
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Fig. 6. Filling-in of Ca II line at 866 nm due to rotational-Raman scattering as a function of solar 
zenith angle for 2 view zenith angles (VZA) typical of minimum (nadir, 0°) and maximum (40°) 
for SCIAMACHY for a relative azimuth angle (AZ) of 90°. 
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Fig. 8. Derived monthly averages for July (left panels) and December (right panels) 2009; Top: 
fluorescence scaled by cos(SZA) (scaled-F) (unitless) from GOSAT (0.696*758 nm + 770 nm); 
Middle: SCIAMACHY 866 nm; Bottom: Aqua MODIS enhanced vegetation index (EVI) 
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Fig. 9. Monthly mean scaled-F derived from 8+ years of SCIAMACHY data (Jan. 2003 - 
July 2011) gridded at 0.5° x0.5° resolution (left) and similarly gridded EVI from Aqua MODIS 
for January through June (top to bottom). 
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Fig. 10. Similar to Fig. 9 but for July through December derived from 8+ years of SCIAMACHY 
data (Jan. 2003 - July 2011). 
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Fig. 11. Regions used for calculating seasonal cycles shown in Fig. 12. 
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Fig. 12. Seasonal cycles derived for the regions shown in the Fig. 11 map, using 0.5° x0.5° 
gridbox monthly means within the region and the vegetation type as listed: Black solid line: 
scaled-F derived from SCIAMACHY (left axis, unitless) for the regions indicated; Horizontal 
lines indicate error of the monthly mean; Vertical lines indicate the standard deviation; Green 
dotted line: MODIS Enhanced Vegetation Index (EVI, right axis, unitless). The numbers indi- 
cated are the number of monthly mean gridboxes used in each plot. 
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Fig. 13. Monthly mean gridded reflectivity near 866 nm as derived from 8+ years of SCIA- 
MACHY data (Jan. 2003 - July 2011) similar to Fig. 9. 
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